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ABSTRACT 

It has recently been proposed, on the basis of eclipse-timing data, that the eclipsing 
polar cataclysmic variable HU Aquarii is host to at least two giant planets. However, 
that result has been called into question based upon the dynamical stability of the 
proposed planets. In this work, we present a detailed re-analysis of all eclipse timing 
data available for the HU Aquarii system, making use of standard techniques used to 
fit orbits to radial- velocity data. We find that the eclipse timings can be used to obtain 
a two-planet solution that does not require the presence of additional bodies within 
the system. We then perform a highly detailed dynamical analysis of the proposed 
planetary system. We show that the improved orbital parameters we have derived 
correspond to planets that are dynamically unstable on unfeasibly short timescales 
(of order 10"^ years or less). Given these results, we discuss briefly how the observed 
signal might in fact be the result of the intrinsic properties of the eclipsing polar, 
rather than being evidence of dynamically improbable planets. Taken in concert, our 
results highlight the need for caution in interpreting such timing variations as being 
planetary in nature. 

Key viTords: 
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1 INTRODUCTION 

Cataclysmic variables (CVs) are interacting binary stars 
composed of a white dwarf primary and a Roche lobe fill- 
ing M dwarf secondary. In the case of HU Aqr, an AM Her 
class CV, the material being accreted by the primary from 
the secondary is channelled along an accretion stream by 
the white dwarf's magnetic field. A comprehensi ve overview 
of these systems can be found in iHellieil (|200l'). Such sys- 
tems are known to exhibit quasi-periodic variations in their 
photometry. A number of factors can contribute to those 
variations, with well accepted causes including the level of 
activity of the secondary (star spots and associated stellar 
cycles) and also that star's shape. 

A number of recent studies have suggested that cer- 
tain CVs host planetary-mass companions. As the postu- 
lated planetary companions orbit the central stars, they 
cause those stars to move back and forth as they orbit 
around the system's centre of mass. As a result, the dis- 
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tance between the Earth and the host stars varies as a func- 
tion of time, meaning that the light from the stars must 
sometimes travel further to reach us than at other times. 
This effect results in measurable variations in the timing 
of mutual eclipse events between the two stars that can 
be measured from the Earth. Using this method, planetary 
mass compan ions have recently been ann ounced around the 
CVs UZ For (iPotter et al.ll201lh. NN Ser (iBeuermann et al ' 
2010), DP Leo l|Qian et al.ll2010[ ). and HU Aqr l|Qian et al. 
I2OIIIV In each of these studies, the eclipse timings are 
first fitted to a linear ephemeris. The residuals from this 
ephemeris (the O — C, or Observed — Calculated timings) 
are then plotted and found to show further, higher-order 
variations. These O — C timings can then be fit with one 
or more superposed Kcplcrian orbits in a manner essentially 
identical to t hat employ ed in Dq ppler radial- velocity planet 
de tection. In Horner et Id] l|201lh . we used the methodology 
of [Marshall et aLn2010l ) to simulate the long-term dynam- 
ical stabi lity of the two g iant planets proposed as orbiting 
HU Aqr lOian et al.ll201ll 'l. We showed that the nominal 2- 
planet solution was extremely unstable on short timescales 
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(~ 10^ yr), unless the outer planet orbited in a retrograde 
and coplanar sense relative to the inner (i.e. with the two 
planetary orbits inclined by 180 degrees to one another). 
Given that such a configuration seems highly unlikely, we 
suggested that either the system is currently undergoing a 
dynamical rearrangement (also highly improbable given the 
^ 10® yr age of the post-main-sequence primary) , or that 
the system is sig nificantly different from that proposed by 
iQian et all (|201l] ). 

In this work, we apply the standard statistical methods 
used by the radial- velocity planet s earch community to the 
timing data for HU Aqr given bv lOian et al.l ()201lh . Sec- 
tion 2 briefly describes the observational data used for our 
analysis. In Section 3, we detail the analysis methods applied 
to these data, and the resulting planetary system configura- 
tions implied. In Section 4 we discuss the dynamical impli- 
cations of our results, before exploring possible alternatives 
to the planet hypothesis in Section 5. 
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2 OBSERVATIONAL DATA 



Following the discovery paper by iQian et all HoTl), we 
make use of the same 82 O — C eclipse egr ess times, con- 
sistin g of 72 data points from the literat ure (ISchwarz et al] 
|2009| ). and 10 new timings presented by iQian et al.l (|201lh . 
All eclipse timings have been fitted with a linear ephemeris 
given by 



2449102.920257 + 0.0868294936B. 



(1) 



The residuals from the fit to this linear ephemeris are plot- 
ted in Figure[T] The root-mean-square (RMS) scatter in the 
timing data is 13.4 seconds, and the reduced of the linear 
ephemeris is quite large at 108.6. There are significant devia- 
tions from the linear ephemeris, suggestive of additional per- 
turbing bodies which result in periodic eclipse timing varia- 
tions. We therefore proceed with fitting Keplerian orbits to 
these signals. 



3 ORBIT FITTING 

In this section, we detail the orbit-fitting proce ss for two 
cases. First, in §3.1, we follow lan et all (|201ll ). and con- 
sider the data after the removal of a quadratic trend. In §3.2, 
we attempt to fit the data without removing a quadratic 
trend. It is well accepted that in fitting radial-velocity data 
(to which results on eclipse timing are clearly analogous) , if 
a long-period object is suspected, the removal of a quadratic 
trend from the data is not ideal. The physically meaningful 
function to fit and remove is a Keplerian (if the variation 
is thought to be due to planetary orbits). As the Keplerian 
function is a complex one, and not necessarily well approxi- 
mated by a quadratic, it is both preferable and more rigorous 
to attempt to fit a Keplerian orbit - even if the parameters 
so derived are not well-constrained, one will at least not 
introduce spurious signals due to a poor match between a 
quadratic and a Keplerian. 



Figure 1. O — C eclipse timings for HU Aqr, after fitting a lin- 
ear ephemeris (given by Equation 1). A quadratic fit to these 
data is ove r plotte d as a solid line. The data are those used in 
iQian et all l201lD. which incl ude 72 observations initially pre- 
sented bv lschwarz et aL 1 l l2009h . At least one sinusoidal variation 
is evident, suggesting the presence of at least one additional per- 
turbing body. 



3.1 Removing a quadratic trend 

In order to test the results of iQian et al.l (|201ll ), we match 
their approach by fitting Keplerian orbits after removing a 
long-term quadratic trend. First, we fit a quadratic trend to 
the original timing data; the fitted parameters are 

(O - C) = (-4.1 X 10"'^)( JD)^ + 0.43(JD) - 1.10 x 10^ (2) 

where JD is the observation date in the form (Julian Date- 
2400000). Then, we fit a single planet (Model Al). A stan- 
dard Lomb-Scargle periodogram [Lomb 1976 ; Scarg k 198^) 
shows a clear signal near 3500 days (Figure [H). We fit a 
Keplerian orbit model using GaussFit ( Jeffervs et al.|[r987l ). 
The single-planet solution is given in Table[T]as "Model Al." 
The residuals to the 1-planet fit are shown in Figure (3] as 
is the periodogram of those residuals. After removing the 
dominant periodicity, there is a significant peak at ^ 8000 
da ys. Using the b ootstrap randomisation method described 
bv lKiirster et al.] ([l997), this peak is found to have a false- 
alarm probability (FAP) of <0.01% (10000 bootstraps). Ow- 
ing to the high significance of this residual peak, and the 
large RMS of the 1-planet residuals (6.57 sec: Table [TJ, we 
proceed by fitting a second planet. 

Given the substantial uncertainty in the system param- 
eters, we used a genetic algorithm to explore a wide param- 
eter space (e.g. Cochran et al. 2007, Tinney et al. 2011). 
The initial range of orbital periods supported by these data 
were first estimated by the periodogram analysis described 
above. The parameters of the best 2-planet solution obtained 
by the genetic algorithm were then used as initial inputs for 
the GaussFit least-squares fitting procedure used above. The 
two-planet fit and a periodogram of its residuals are shown 
in Figure [5] The parameters of the 2-planet fit are given in 
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Figure 3. Results of single-planet fit on data with quadratic trend removed (Model Af). Left panel: Data and model fit for a single 
planet. The quadratic trend we have removed is superposed on the f-planet Keplerian model (solid line). Right panel: Periodogram of 
residuals after fitting for and removing the dominant signal, a Keplerian orbit with period 3538 days. A significant peak is seen at a very 
long period. 




Figure 2. Periodogram of the O-C timing data for HU Aqr, using 
raw data (residuals from a linear ephemeris) with quadratic trend 
also removed. A strong signal is evident near 3500 days. 

Table [1] as "Model A2." Since the total duration of the data 
set is 6118 days, and the best-fit period for an outer body 
is 7215 days, there remains significant uncertainty in the 2- 
planet fit. In particular, the Keplerian orbit-fitting process 
failed to converge when the outer planet's eccentricity and 
periastron argument (cj) were allowed to be free parameters. 
Hence, the values for these parameters shown in Table[T]were 
held fixed at the best result from the genetic algorithm. As 
shown in Figure [G] these two parameters are almost com- 
pletely unconstrained. This is not surprising given that the 



Figure 4. Periodogram of the O-C timing data for HU Aqr, 
using the raw data (residuals from a linear ephemeris) with no 
additional quadratic trend removed. A strong signal is evident 
near 5000 days. 

period of the outer planet is larger than the entire length of 
the data set. 

3.2 No quadratic trend 

In this subsection, we explore the possibility that the re- 
moval of a quadratic trend has confounded the orbit-fitting 
process by absorbing signal due to a long-period outer 
planet. Here we repeat the fitting procedures as above, but 
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Figure 5. Left panel: Two-planet fit (Model A2) on data with quadratic trend removed. The quadratic trend we have removed is 
superposed on the 2-planet Keplerian model (solid line). Right panel: Periodogram of residuals to this fit. No further significant periods 
are evident. 
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Figure 6. Results of genetic algorithm fitting for 2 planets (Model A2). Left panel: Dependence of reduced o^i the outer planet's 
eccentricity. Right panel: Same, but for the outer planet's argument of periastron (oj). These parameters are essentially unconstrained, 
as nearly the entire allowed range is within 1.0 of the minimum. 



using the HU Aqr data which have not had a quadratic 
trend subtracted. First, we considered a single planet by 
performing a periodogram search (Figure [4]), which shows 
the dominant periodicity to be at 5000 days. However, the 
standard approach of fitting a Keplerian orbit with GaussFit 
failed. As the reduced of the best-fit genetic algorithm 
result was an inordinately high 21.2, we attribute the fail- 
ure of the least-squares method to the presence of additional 
signals. In Table [U we give the best-fit results for a single 
planet ("Model Bl") from 100,000 iterations of the genetic 
algorithm. One-sigma uncertainties are estimated by noting 
the change in each parameter required to increase the re- 



duced by 1- The dependence of x^ on each parameter is 
shown in Figure [T] As with the one-planet solution in the 
previous trial (data which included a quadratic trend), the 
RMS scatter about a 1-planet model is quite large at 8.0 
seconds, compared to the mean measurement uncertainty of 
2.4 seconds. The fit and a periodogram of its residuals are 
shown in Figure [HI it is clear that one planet is not sufficient 
here. The highest periodogram peak is at a period of 2128 
days, again with a FAP<0.01%. We thus proceed to fit a 
second Keplerian orbit. 

Again we employ the genetic algorithm to explore the 
wide and uncertain parameter space for a 2-planet model 
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Figure 7. Results of single-planet genetic algorithm fit for HU Aqr (Model Bl). Each panel shows the dependence of reduced on a 
particular planetary parameter. The uncertainty of each parameter is estimated by the range over which the reduced increases by 1 
from the minimum. 



("Model B2"). First we examine the short-period option, 
as prompted by the periodogram results in Figure [8l We 
ran the genetic algorithm for 100,000 iterations, then at- 
tempted a standard least-squares fit on the best result. As 
with Model Bl, no such fit would converge with the ec- 
centricities and periastron arguments of the two planets 
free. The best-fit model with a short period for the sec- 
ond planet (actually making it the inner of the two planets) 
resulted in a reduced of 4.06 and an RMS of 5.97 sec- 
onds. The planetary parameters r esulting from this fit are 
quite similar to those proposed bv lQian et all (|201ll ). with 
Pinner = 1947±10 days and Pouter = 4429±113 days. How- 
ever, this is substantially worse than the two-planet fit from 
Model A2, and also worse than the long-period fit obtained 
below. 

Allowing the genetic algorithm to choose long periods 
for the second planet, we obtain a far better solution, with 



a reduced of 0.80. The parameters of this fit are given 
in Table [T] as Model B2; this fit and a periodogram of its 
residuals are shown in Figure [9] Both Models A2 and B2 
support a long period for the second planet, so the short 
period case discussed briefly above is rejected. Again, the 
fitting process failed when e and lu for the outer planet were 
free parameters, so we fixed their values at the best-fit from 
the genetic run (100,000 iterations). Uncertainty estimates 
are obtained from the plots in Figure [TOl for uj, the surface 
has two minima, and so no formal uncertainty is quoted - 
this parameter is very poorly constrained due to the long 
period of the outer planet. 



In iQian et al] (|201ll ) and iHorner et al.l (|201ll) . it was 

suggested that a third, distant outer planet may be present 
in the HU Aqr system. However, in light of the results of the 
2-planet fits given in this section, which feature reduced 
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Figure 8. Results of single-planet fit for HU Aqr (Model Bl). Left panel: Data and model fit for a single planet. Clearly a single planet 
is inadequate to fit these data. Riglit panel: Periodogram of residuals after fitting for and removing the dominant signal, a Keplerian 
orbit with period 4728 days. A significant peak is seen near 2100 days. 
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Figure 9. Left panel: Two-planet fit (Model B2) to the observed data with no quadratic trend removed. Right panel: Periodogram of 
residuals to this fit. No further significant periods are evident. 



less than 1.0, we see no need to invoke additional bodies to 
adequately fit the available data. 



4 DYNAMICAL ANALYSIS 



In summary, our analysis of two slightly different ver- 
sions of the HU Aquarii data yields evidence for two plan- 
ets: a moderately well constrained one at P = 4647 — 4688d 
with e ~0.2, and a somewhat more poorly constrained one 
at P = 7215 — 8377d with a poorly constrained but non- 
zero eccentricity. While the outer planet's period varies by 
~ 1200d between these two solutions, we note that this rep- 
resents just a 2-a difference, given the period uncertainties. 
Hence, there is a long period outer signal present, even if its 
period is not well determined. 



Following the results detailed above, the observational data 
yield two distinct 2-planet solutions which are essentially 
identical in terms of their goodness-of-fit criteria. In both 
solutions, the best-fit syste m parameter s are significantly 
different from those given bv lQian et al.l l|201ll), whose dy- 
namic al stability was studied in some detail in Horner et al.l 
(j20i3)- 

In that work, the authors argued that the extreme 
levels of dynamical instability displayed by the planetary 
system provided firm evidence that, at the very least, the 
true parameters of the planets in the system were greatly dif- 
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Figure 10. Results of genetic algorithm fitting for 2 planets (Model B2). Left panel: Dependence of reduced on the outer planet's 
eccentricity. Right panel: Same, but for the outer planet's argument of periastron (oj). This parameter displays two minima, and so is 
very poorly constrained by the available data. 



ferent to those obtained through the analysis of lOian et al.l 
l|201ll ). 

How do the new solutions for the proposed HU Aqr 
planets stand up to the same test? In order to closely exam- 
ine the dyn amical stability of th e ne wly proposed solutions , 
we followed [Horner et al.l l|201ll 'l and lMarshall et al.l (|2010l ') 
and performed highly detailed dynamical simulations of the 
potential planetary systems. Such simulations serve as a crit- 
ically important additional test, since the algorithms used 
to this point include no physics; rather, the Keplerian fitting 
methods are simply seeking a lowest-x^ solution regardless 
of the phys icality of the r esultin g system parameters. 

As in iHorner et al.l (|201ll '). we used the Hybrid in- 
tegrator within the A'^-body dynamics package Mercury 
IChambers Ill999l ) to perform our integrations. For the two 
scenarios in question (Model A2 and Model B2), we cre- 
ated 50625 test planetary systems. In each case, we followed 
our earlier work, and kept the initial orbit of the innermost 
planet fixed at its nominal best fit value. The initial orbit of 
the outermost planet was then varied systematically in semi- 
major axis a, eccentricity e and mean anomaly M, such that 
a total of 50625 unique solutions were tested. For our tests 
of each of the two models, 45 initial values of a were tested, 
spread evenly across the full ±3(t error range in that pa- 
rameter. Similarly, 45 initial values of e were tested in each 
case, with 25 different AI being considered for each a — e pair- 
ing. For Model A2, in which the e of the outermost planet 
was unconstrained, we tested eccentricities ranging between 
0.005 and 0.995, whilst for Model B2, the tested e values 
were sprea d across the full ±3 a errors given in Table [1] 

As in iHorner et all l|201ll '). we followed the dynamical 
evolution of each test system for a period of 100 million 
years, and recorded the times at which either of the planets 
was removed from the system. Planets were removed if they 
collided with one another, hit the central body, or reached 
a barycentric distance of 100 AU. 

The results of our dynamical integrations can be seen 



in Figure [TT] and Figure 1121 The 3 — a region around the 
nominal orbits in a — e space is clearly highly dynamically 
unstable. As was seen in the dynami c al stu dy of the plan- 
etary system proposed bv lQian et al.l l|201ll ). low eccentric- 
ity orbital solutions which keep the planets separated at all 
times by at least three Hill radii offer some moderate in- 
crease in the potential stability of the system. We note that 
these regions lie far from the central 1 — cr of the error el- 
lipse. Furthermore, a fit with the outer planet's period and 
eccentricity fixed in this stable region has a reduced of 
14.9, far worse than the best-fit results given in Table [1] 

Muc h as was the case for the planetary system as pro- 
posed by iQian et al.l (|201ll ). we therefore find that both of 
the best-fit models detailed above (Model A2 and Model B2) 
fail to stand up to dynamical scrutiny. We find that the re- 
sulting dynamical instabilities make it exceedingly unlikely 
that the eclipse-timing variations observed for HU Aqr are 
truly the result of perturbations from planetary mass ob- 
jects in the system. It seems reasonable, therefore, to ex- 
amine more closely whether any other effects could cause 
timing variations of an appropriate periodicity and scale in 
such eclipsing polar systems. 



5 ALTERNATIVES TO THE PLANET 
HYPOTHESIS 

There are many causes of variability, spanning timescales 
from seconds to years, which are inherent to the nature of 
CVs. All of these mechanisms may have an observable im- 
pact on the variation of the O ~ C curve, which may act to 
dilute, obscure, or mimic a signal that would otherwise be 
attributable to the presence of exoplanet(s). 

At the time scale considered here (i.e. thousands of 
days) the most likely source of the observed variation in the 
O — C curve is the secondary star in the system, the M dwarf. 
With a rotation period on the order of a few hours (the result 
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Figure 11. Left panel: Median lifetime of the outermost planet in the HU Aqr system, based on the orbits fit by Model A2, as a 
function of the planet's semi-major axis a, and eccentricity e. The lifetime shown in each square is the median of the 25 distinct runs 
(each of which had the planet start at a different initial mean anomaly, Af). In both panels, the vertical dashed lines indicate relevant 
mean-motion resonances. This panel reveals that the whole phase-space of potential orbits described in Model A2 is highly dynamically 
unstable, with median lifetimes typically much less than 10® years. Right panel: Same as the left-hand panel, except that the mean, 
rather than the median, of the 25 initial mean anomaly values is calculated for each bin across the plot. Whereas the left panel clearly 
showed the overall lack of stability of the plausible orbital solutions given by Model A2, plotting the mean lifetimes of each bin reveals 
that, for a small region of orbital space for which the outermost planet could approach no closer to the innermost planet than 3Rh, 
some fraction of the tested orbits were dynamically stable on much longer timescales. 
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Figure 12. Left panel: Median lifetime of the outermost planet in the HU Aqr system, based on the orbits fit by Model B2, as a 
function of the planet's semi-major axis a, and eccentricity e, as described above. In both panels, the vertical dashed lines indicate 
relevant mean-motion resonances. The great bulk of the tested a — e phase space is highly dynamically unstable. Right panel: Same as 
the left-hand panel, except that the mean, rather than the median, of the 25 initial mean anomaly values is calculated for each bin across 
the plot. As was the case for the results of Model A2, this reveals that a small fraction of orbits tested can display long-term dynamical 
stability, though such orbits remain in the minority. 



of it being tidally locked in its rotation about the primary, 
the white dwarf), the dynamo effect in the secondary could 
be large. Assuming that the stars in such systems have a 
magnetic cycle similar to the Sun's double-peaked 22-year 
cycle, then their angular momentum distribution will change 
over time. This will have the effect of altering the shape of 
the secondary, which in turn affects the gravitational at- 



traction between t he primary and secondary - and hence 
the orbital period (|Warneij|l98i lApplegatdllOQ^ ) . 

This effect has already been observed in several CVs. 
Examples include U Gem, which displays a ~1 minute 
variation in orbit al period over an eight year time scale 
(lEason et al ]| 19831 ). EX Dra, whose period varies by 1 . 2 min - 
utes over approximately four years l|Baptista et al.l |2000| ) , 
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Table 1. Orbital Solutions for HU Aquarii 



Model 


Parameter 


Inner Planet 


Outer Planet 


Al 


Orbital Period (days) 


3538±54 






Amplitude (sec) 


12.5±0.8 






Eccentricity 


0.26±0.05 






uj (degrees) 


171±12 






To (JD-2400000) 


53092±517 






Urbitai Haams (AU ) 


4.DD±U.li 






Msini (Mjup) 


6.08±0.13 






xl 


7.40 






RMS (sec) 


6.57 




A2 


Orbital Period (days) 


4647±36 


7215±603 




Amplitude (sec) 


20.2±0.5 


23.4±1.9 




Eccentricity 


0.19±0.02 


0.53 (fixed) 




UJ (degrees) 


166±5 


190 (fixed) 




To (JD-2400000) 


53074±50 


58060±599 




Orbital Radius (AU) 


5.59±0.11 


7.50±0.50 




Msini (Mjup) 


8.19±0.24 


7.07±0.41 




xl 


0.69 






RMS (sec) 


2.50 





Bl Orbital Period (days) 

Amplitude (sec) 
Eccentricity 
LJ (degrees) 
To (JD-2400000) 
Orbital Radius (AU) 
M sini (Mjup) 



4728l^»» 
15.1±1.6 
0.09+0"5 

181 



-10 

53991±60 
5.66±0.31 
6.05±0.27 



B2 Orbital Period (days) 4688±177 

Amplitude (sec) 14.0±2.1 

Eccentricity 0.20±0.04 

UJ (degrees) 254±12 

To (JD-2400000) 53640±116 

Orbital Radius (AU) 5.62±0.22 

Msini (Mjup) 5.65±0.20 



8377±610 

27.9±3.7 

0.38+55;^^ (fixed) 

332 (fixed) 

60126±612 

8.28±0.50 

7.64±0.12 



xl 



RMS (sec) 



0.80 
2.49 



of planet formation, accreting from some of the mate- 
rial shed by the primary star ijTavani fc Brookshawl Il992l : 
|Phinney fc Hansen|[l993l : [Hansen et al.ll2009l )." This "second- 
generation" planet formation scenario was explored after 
the discovery of the terrestr ial-mass planets orbiting pul- 
sars (jWolszczan fc Fr"ailll992l ). However, the mechanism by 
which Jupiter-mass planets could form in such an evolved 
system is not currently clear. 



6 SUMMARY AND CONCLUSIONS 



xl 


21.2 


RMS (sec) 


8.00 



and EX Hya, whose perio d is modulated on a tim escale of 

approximately 17.5 years (|Hellier fc Sproat3ll992| ). 

On the question of planetary survivability, Qian et al] 

mentioned that a gas giant planet could survive 
the planetary nebula stage, so long as it was located be- 
yond orbital distances of a bout 3 AU l|Villaver fc Livij 
l2007l : iKunitomo e^al]|2011^ . However, to make it to the 
planetary nebula phase, the planet must first survive the 
Asymptotic Giant Branch (AGE) phase. A normal AGB 
star with a main sequence mass of 5 Mq (extrapolated 
from the white dwarf mass of 0.88 Mq) is expected to 
reach a radius of 5. 25 AU, envelopin g both of the pos- 
tulated planets from lOian et al.l (|201ll ). A planet entering 
the envelope of an AGB star would clearly have very lit- 
tle chance of surviving. Planets orbiting post-main-sequence 
stars are thought to be either survivors of the planetary neb- 
ula o r supernova phase ( Colgateiri970l : iPostnov fc Prokhorovl 
1 19921 : IVeras et all bOlll T or formed from a second phase 



We have revisited the work of iQian et al. I ('2011'), who re- 
ported that timing variations in the mutual eclipses between 
the primary and secondary components of the HU Aqr sys- 
tem were evidence that there were at least two Jupiter-mass 
planetary companions orbiting around the two central stars. 

In this work, we applied the key tools used in the de- 
tection of planets around main sequence stars t hrough the 
radial velocity technique to the data used by iQian et al.l 
(l201ll ). Our analysis resulted in two distinct two-planet fits 
which could, in theory, explain the observed timing varia- 
tions. 

Our derived orbital solutions differ significantly from 
those presented in that earlier work, but still fall prey to the 
same dynamical drawback. Simply put, the planets neces- 
sary in order to explain the eclipse-timing variations prove 
dynamically unstable on timescales so short as to seem un- 
feasible. 

Our results therefore suggest that some other mecha- 
nism must instead be invoked in order to explain the ob- 
served variations. The most likely candidate, based on ear- 
lier studies of eclipsing polar systems, is that the observed 
variations are the result of the interaction between the mag- 
netic fields of the stars in the course of their stellar cycles. 
Such variations have been observed in a number of eclipsing 
polar systems in the past, and would be expected to yield 
variations on the scale observed, over similar timescales. 

In light of these results, it would seem prudent in fu- 
ture to consider such behaviour as a potential source of sig- 
nals that could mimic the presence of planets orbiting cat- 
aclysmic variable stars. Rigorous dynamical testing of any 
planetary systems resulting from the analysis of transit tim- 
ing data for such stars should become a key component of 
the analytical process. These tests will prove critical in dis- 
tinguishing between CVs which might plausibly host such 
interesting planetary systems, and those in which such plan- 
ets are all but impossible. 
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Table 2. Eclipse Timing Residuals for HU Aqr 



JD- 2400000 


O-C (s) 


Uncf 


49217.436369 


-22.2 


1.0 


49217.523190 


-23.1 


1.0 


49217.610010 


-22.6 


1.0 


49217.696830 


-19.8 


2.0 


49218.651855 


-21.2 


2.0 


49218.738675 


-20.4 


2.0 


49221.603749 


-21.4 


2.0 


49221.690569 


-21.6 


2.0 


49221.777389 


-20.9 


2.0 


50325.959283 


-2.5 


2.4 


50326.046103 


-2.5 


2.4 


50338.895523 


-3.7 


2.4 


50285.414154 


-4.3 


10.0 


50285.500975 


-8.3 


10.0 


50286.455999 


-9.4 


10.0 


50328.390254 


6.0 


10.0 


50328.477074 


7.6 


10.0 


50330.387123 


-7.1 


10.0 


50330.473944 


-6.7 


10.0 


51081.383615 


8.3 


10.0 


51481.278394 


20.9 


2.0 


51703.625447 


22.3 


7.0 


51704.580472 


20.1 


7.0 


52145.367661 


22.8 


4.0 


51703.625447 


22.3 


8.0 


51704.580472 


20.1 


8.0 


51821.440735 


24.7 


1.0 


51821.527555 


24.7 


1.0 


51350.874147 


15.1 


1.0 


51353.826041 


15.4 


1.0 


51354.867886 


16.6 


1.0 


51731.494797 


23.5 


1.0 


51731.581617 


23.7 


1.0 


52174.278855 


20.9 


1.0 


52174.365676 


21.1 


1.0 


52411.559018 


14.9 


1.0 


52552.381713 


11.3 


1.0 


52553.336737 


8.9 


1.0 


52787.665007 


3.0 


1.0 


52789.575056 


1.7 


1.0 


53205.444789 


-6.8 


1.5 


53205.531609 


-7.1 


1.5 


53209.525348 


-6.0 


1.5 


53212.564062 


-5.7 


1.5 


53293.307038 


-7.6 


1.5 


53295.303907 


-9.2 


1.5 


53296.258931 


-7.5 


1.5 


53299.297645 


-7.6 


1.5 


53533.539094 


-6.8 


1.5 


53918.500764 


-4.8 


1.5 


53925.446396 


-4.9 


1.5 


54320.479233 


1.5 


1.5 


54320.566053 


-0.3 


1.5 


52411.211737 


11.6 


5.0 


52779.937991 


3.7 


5.0 


52411.385378 


16.8 


2.0 


52411.385378 


17.3 


2.0 


52411.472198 


16.8 


2.0 


52413.642708 


12.3 


1.0 


52414.684553 


16.2 


1.0 


52783.671268 


2.6 


1.0 


53504.888361 


-5.8 


0.5 


53505.843386 


-5.9 


0.5 


53505.930206 


-5.8 


0.5 


53506.798410 


-5.9 


0.5 



